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Abstract 

To test a multicluster approach for halo nuclei, we give a unified description for 
the ground states of ^He and ^He in a model comprising an a cluster and single- 
neutron clusters. The intercluster wave function is taken a superposition of terms 
belonging to different arrangements, each defined by a set of Jacobi coordinates. 
Each term is then a superposition of products of gaussian functions of the individual 
Jacobi coordinates with different widths, projected to angular momenta Z = or 1. 
To avoid excessively large dimensions and "overcompleteness" , stochastic methods 
were tested for selecting the gaussians spanning the basis. For ^He, we were able to 
calculate ground-state energies that are virtully exact within the subspace defined by 
the arrangements and I values, and we found that preselected random sets of bases 
(with or without simulated annealing) yield excellent numerical convergence to this 
"exact" value with thoroughly truncated bases. For *He good energy convergence 
was achieved in a state space comprising three arrangements with all / = 0, and 
there are indications showing that the contributions of other subspaces are likely 
to be small. The ^He and ^He energies are reproduced by the same effective force 
very well, and the matter radii obtained are similar to those of other sophisticated 
calculations. 



* Permanent address. 



1 Introduction 



The discovery [|I| of thick neutron clouds or "halos" in the surface of some hght nuclei 
has led to a revival of interest in these nuclei, as is shown by numerous experimental 
and theoretical studies (see e.g. |^). All candidates for halo structure are in the p-shell 
region, where systematics can be studied with the shell model, but the shell model runs 
into difficulties just in producing halos. To bring about a halo with a simple-minded en- 
largement of an harmonic-oscillator (h.o.) basis would in most cases require prohibitively 
large dimensions. If, alternatively, the single-particle orbits are chosen so as to exhibit the 
loosely-bound character of the halo [Q, then the treatment of the centre-of-mass (cm.) 
motion, of the Pauli principle and of continuum orbits pose extra problems. 

The cluster-orbital shell model (COSM) |^ has solved these problems successfully, and, 
at least qualitatively, reproduced some known properties of the halo nuclei ^He 0, ^, ®He 
1^ and ^^Li But the COSM is not necessarily the ultimate model either. Even for the 
simplest case of ^He, it is difficult to achieve convergence 0. With the number of valence 
particles increased, this problem gets aggravated. On the other hand, the nuclei that have 
two halo nucleons are often "borromean": removal of one nucleon prompts a breakup of 
the rest by emission of another nucleon. It was shown that the three-body dynamics 
characteristic of such a system can be better described by inclusion of a core+n+n type 
three-cluster component 0. 

The alternative models that focus just on the three-body aspects are the three-body 
models (see e.g. P, ||). None of the three-body calculations has, however, achieved 
unquestionable success. One problem is that the core-nucleon interaction, at least its 
off-shell behaviour, is not well-known. The accuracy of the treatment of antisymmetriza- 
tion through Pauli projection (let alone through inclusion of a repulsive interaction |llO| ) 
depends crucially on the off-shell behaviour of the potential. Even if this treatment is 
satisfactory for the solution of the dynamical problem, there are extra three-cluster Pauli 
effects in the spectroscopic amplitudes characterizing the disintegration [^. Moreover, 
it would be illusive to believe that cores like ^Li are perfectly inert. 

All these shortcomings call for fully microscopic treatments that preserve some el- 
ements of the few-body approaches as well. Although microscopic cluster models are 
general enough and their techniques are highly sophisticated |TB|, they have not 
been used too much to describe systems containing more than two clusters. It has been 
demonstrated recently for an a+p+n model of ^Li that a microscopic approach, in which 
the treatment of the relative motion is comparable with the three-body approaches is fea- 



sible ||15|, and similar calculations for the He- Li- Be isospin triplet have been published 



in another paper [0]. 

In the present paper we shall apply a multicluster approach to the ground state (g.s.) 
of the three-cluster system ^He=Q;+n+n and of the five-cluster system ^He=Q;+n+n+n+n. 



The present approach is similar to the previous one |T^, except for some practical sim- 
plifications, which are compelling for ^He, and, for consistency, are adopted for ^He as 
well. The treatment of ^He is a prerequisite of that of ^He because some parameters are 
to be fitted to ^He. The case of ^He is used, at the same time, as a testing ground of 
the approximations introduced. We have solved the dynamical problems of these nuclei, 
hence obtaining the g.s. energies, and, to test the wave functions, we have calculated the 
rms point matter radii. Further applicaltions of our wave functions will be published in 



2 



a forthcoming paper. 

The model is based on a hnear variational method. The wave function is a superposi- 
tion of cluster-model terms belonging to different sets of Jacobi coordinates or "arrange- 
ments", and the only parameters to be varied are the combination coefficients. Unlike 
before [0, we now choose purely central effective nucleon-nucleon interactions. This 



leads to a less realistic description of the relative motion, but it is more consistent with 



a purely s-state model of the a intrinsic motion |15]. In consistency with the use of a 



central force, we allow for fewer values of the intercluster orbital angular momenta. This 
is expected to be a fair approximation because the motion of two halo neutrons relative 
to each other, as well as of neutron pairs relative to a core, is dominated by s waves. 
Just as before, we expand the relative-motion functions in terms of gaussians of different 
widths. To counter the explosion of the basis size, we now have to introduce special trun- 
cation schemes for these expansions. We have examined systematic as well as stochastic 
methods for constructing the basis. As will be seen, the stochastic methods have proved 
more successful. It is owing to the stochastic variational procedure that this extremely 
versatile model is feasible for practical purposes. 

The plan of the paper is as follows. In sect. 2 we outline the model used. In sect. 3 
we introduce various basis truncation schemes and test them for ^He. Sect. 4 comprises 
the results for ^He, and sect. 5 contains a summary and some conclusions. 



2 The multicluster model 

The model we use is a generalization of the variational multiconfiguration multicluster 



approach outlined in ref. [T^, and involves an a. cluster and n single-neutron clusters. It 
employs a trial function, which is a sum over various cluster arrangements /x, each asso- 
ciated with a particular set of intercluster Jacobi coordinates {p^*, P^}- For a-fn-t-n, 
there are two non-trivially different arrangements, viz. yU=(a;n)n and /i'=a;(nn), and the 
corresponding intercluster Jacobi coordinates are 

Pi = ri-rc,, P2 = r2 - r(4r„ + ri). 



5' 

p; = r2-ri, P2 = |(ri + ra) - r^,, (1) 



respectively, where r^^ is the position vector of the a cm. and rj are those of the extra 
neutrons. As is seen in fig. |l]a, the Jacobi coordinates for (Q;n)n and a(nn) form a Y- 
shaped and a T-shaped pattern, respectively, and they will be referred to as such. For 
the system a+n+n+n there are four arrangements, [(Q;n)n]n, [a;(nn)]n, (an)(nn) and 
a[(nn)n], while for a+n+n+n+n there are nine: {[(Q;n)n]n}n, {[a;(nn)]n}n, [(an)(nn)]n, 
{a[(nn)n]}n, [(an)n](nn), [Q;(nn)](nn), (an)[(nn)n], a{[(nn)n]n} and a[(nn)(nn)]. Fig. |l|b 
shows a few of the five-cluster arrangements. As an illustration, here we give the Jacobi 
coordinates belonging to the TY-shaped arrangement {[a(nn)]n}n: 

Pi = r2-ri, P2 = K^i + ^2) - Ta, 

P3 = - |(4ra + ri + r2), p4 = r4 - i(4ra + ri + ra + rg). (2) 

The nucleon spins are coupled to 5, and the orbital angular momenta (i = 1, n), 
belonging to the Jacobi coordinates Pj, are coupled to L. A term pertaining to arrange- 
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ment /i and angular momenta [S, JM can be written as 

^'s,il...U)L = ^ { [^5A'^,...„).(P^ P'n)] J J , (3) 

where A is the intercluster antisymmetrizer, ^sms is a vector-coupled product of the a 
intrinsic wave function and the n neutron spin-isospin functions $\ and A')-'^^ ^^-^^ is 
a vector-coupled product of the intercluster relative functions Xikm^iPi) ■ The sequence 
of coupling is chosen so as to follow the pattern of the Jacobi coordinates. The intrinsic 
wave function of the g.s. of the a particle, is an h.o. Slater determinant divided 
by the Os h.o. wave function describing the zero-point vibration of the a cm. The h.o. 
size parameter u = muj/2h is determined by minimization of the a intrinsic hamiltonian. 
(The spatial factor of such a function is equal to a product of Os h.o. functions of three 
intrinsic Jacobi coordinates.) 

Representing each relative motion belonging to any single set of Jacobi coordinates 
by a complete set of functions, one can get a complete basis. The intercluster functions 
XiimiPi) ^^6' however, approximated by finite sums. 



XilmiPi) = ^ '^ikl^iklmiPi) ^ (4) 
k=l 



where rf^.^^ are h.o. eigenfunctions of radial node number zero: 



^iklmiPi) 



v^(2/ + l)!! 



{pfyeM-'^tkipfy]Yimm- (5) 



The size parameters {z/^'^, k = 1, ...,iVf} should be chosen such that the functions Ff^/^ 
adequately span the state space in the relative motion defined by the radial Jacobi coor- 
dinate . This can always be achieved since any well-behaved functions can be approxi- 
mated with any prescribed precision by a combination of such functions. Because of the 
loosely-bound nature of these nuclei, more than one partial wave k may be needed for 
each relative motion, with a rather large N-^ (~ 10) for the radial motion in each partial 
wave, which entails an enormously large basis. Each set of Jacobi coordinates implies a 
particular asymptotic configuration, which involves few partial waves just when described 
in terms of these particular Jacobi coordinates. Therefore, the dimension may be reduced 
by combining the different arrangements: 

^(5,L)JM = E E %,ih...l„)L]JM^ (6) 
M li...l„ 

where the lack of summations over S and L anticipates the choice of a purely central 
interaction. 

The sets of Jacobi coordinates differ in what physical effects they emphasize. The 
Y-shaped and T-shaped branches of the diagrams of fig. tend to magnify the corre- 
lation of the neutrons with a "core" and with each other, respectively. Thus, for ^He 
the arrangement {[(an)n]n}n is best suited to describing a shell-model-like configuration. 
The arrangements [a(nn)](nn) and a[(nn)(nn)] prefer two dineutron clusters, (Q;n)[(nn)n] 
and {a[(nn)n]}n favour the formation of a three-neutron cluster, a{[(nn)n]n} seems to 
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accommodate a four-neutron cluster, while {[a(nn)]n}n, [(an)(nn)]n and [(an)n](nn) are 
moulded to contain one dineutron and two neutrons moving around a core in a single- 
particle manner. It is a priori obvious that some of these configurations are very similar 
to each other, for instance, the last three arrangements only differ in what the core is. 
Thus it seems sufficient to include a few of these configurations. Even so, the basis sub- 
sets belonging to different arrangements do have large overlaps especially for 
small distances, which not only allows but, to avoid (numerical) overcompleteness, also 
necessitates truncations. 

Both the angular momenta and the expansion of the radial functions are to be trun- 
cated. From among highly overlapping subspaces, we shall keep those with pure k = 
because, for several clusters, the calculation of the matrix elements with non-zero k is 
very lengthy. To introduce a truncation of the expansion of the radial function X, we 
need to substitute into X involved in eq. (^: 

^{h...l„)LA — Z^^K,h...ln Ikih---'- nknln\ij>^^ y ) 

K 

where K = {ki...kn}, D'^i_^ = Cikiii---Cni:„i„y the summation over K may be 
truncated such that it will no longer be equivalent to a sequence of summations over 
ki,...,kn- 

The only variational parameters in the model are the coefficients Since z/^^ are 

not varied, the bases will be somewhat redundant in the sense that there will exist smaller 
bases that would produce the same accuracy. To keep the variational method linear is, 
however, advantageous because the solution is reduced to a diagonalization over a non- 
orthogonal basis, and the method is founded on a genuine minimum principle. This implies 
that any estimate for any bound-state energy gives an upper bound, and the energies are 
lowered monotonously as the basis is extended. Moreover, if the basis is extended so that, 
in a limiting case, it fully spans a well-defined subspace of the state space, the energies 
will converge to the exact energies in that subspace. In the foregoing applications whole 
sets of variational calculations will be performed with stepwise extensions of the basis to 
make estimates for such "exact" energies by producing numerical convergence. 

We chose the central effective interaction introduced by Thompson, LeMere and Tang 



Vq , which has the form 

Vit,j) = [Vniu,) + i(l + Pr.mr,,) + 1(1 - Pr,)Vsir,,)] i [n + (2 - u)Pr 

where rij = \rj — rj|, P^j, P[j are the spin- and space-exchange operators, and V{rij) are 
of gaussian shape, and we neglected the Coulomb force, which would just shift the energy 
scale. This interaction reproduces, at least qualitatively, the most important low-energy 
nucleon-nucleon phase shifts, which implies that it does not bind the dineutron. On the 
other hand, by tuning the exchange parameter u, we can make the a+nucleon system 
comply with the borromean character of ^He. The parameter u can be adjusted so as to 
reproduce the P3/2 a+n phase shift and the ^He binding energy E with respect to the 
three-body breakup threshold simultaneously, with reasonable accuracy. Without having 
a spin-orbit force, this is the best one can do. Since in ^He the P3/2 shell can accommodate 
all valence neutrons, this prescription is expected to cause little error. We chose m=1.14, 
which produces E = —1.016 MeV, while the experimental value is —0.975 MeV. The a 
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size parameter, rms point matter radius and energy implied by this force are v = 0.303 
fm~^, 1.363 fm and = —25.581 MeV, respectively. 

The matrix elements involving the basis functions are expressed as integral transforms 
of matrix elements of Slater determinants, of shifted Gaussians, projected to total angular 



momentum JM. The technique applied is similar to that used in our earlier work |T5 



3 Description of ^He 

In the description of a bound-state three-cluster system, like ^He, the use of a large enough 
basis poses no problem. For the g.s. of ^He, a central force renders the quantum numbers 
to pure S = L = and h = h- Appreciable contributions only come from the subspace 
defined by the T-shaped arrangement with (hjh) = (0,0) ("Too subspace") and by the 
Yoo and Yn subspaces defined analogously. Radial bases of the types of eq. (||) with 
parameters 

«)"'^' = 0.7 X 1.4'=^ fm ih = l,..., 10; z = 1, 2) (9) 

yield energies that are virtually exact within the model space defined by the arrange- 
ment(s) and value(s) of (/i, I2) adopted. In table |l| we give these energies along with the 
corresponding rms point matter radii for the three subspaces Too, Yoo and Yn and their 
unions. 

The energies obtained in the individual subspaces are very close to each other, which 
shows that they are almost equally important, and yet, when they are combined, the 
binding gets stronger moderately. This behaviour is typical of strongly overlapping sub- 



spaces |T^, |I6|. Although, as can be inferred from the table, none of the subspaces 
is large enough by itself, the combination {Too,Yii} does look like that. The fact that 
the energy produced by state space {Too,Yii} is hardly pushed down by inclusion of Yoo 
(with the radius unchanged) indicates that the "full" wave function generated in the space 
{Too,Yoo,Yii} is in fact almost entirely contained in subspace {Too,Yii}. 

The observed overlap between the subspaces suggests that the basis is largely redun- 
dant. It would be obviously desirable to reduce its dimension by pinpointing redundant 
elements and dropping them. This is, however, not a trivial task because the overlaps of 
each particular element with all others are diverse. In the following we shall test various 
tactics of omitting elements and choosing elements so as to produce less redundant bases. 
We should emphasize that, even when the dimension is large (e.g., 300), the bulk of the 
computing time is spent on the computation of the matrix elements rather than on the 
diagonalization. Therefore, unlike in an earlier work on basis optimization [^, our aim is 
to reduce the number of matrix elements to be computed, and to this end, we are willing 
to make more diagonalizations. 

The basis of each subspace defined by eq. (^ can be depicted in a (61, 62) = ((^f)^^^^, 
(z/g)^^/^) sheet as is shown in fig. 0a. We see that in each subspace for each Jacobi 
coordinate each function F occurs 10 times while the u for the other Jacobi coordinate 
runs over all values. It is thus natural to ask whether one could get a satisfactory basis 
by a straightforward "thinning" . Restricting ourselves to the Too space, we discarded half 
of the basis elements in a chessboard-like pattern shown by the open circles in fig. ^ (i.e. 
we kept {F^'fc^j^, ki = 2,4, 10} for k2 odd and {F^^^^^, ki = 1,3, ...,9} for ^2 even). We 
obtained E = —0.312 MeV, which implies too much loss (0.083 MeV) of binding energy. 
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Alternatively, one may set up the basis with a stochastic method. First we present 
calculations with randomly chosen subsets of the basis elements defined by eq. as 
exemplified by fig. ^d. Fig. ^ shows the energy versus dimension curves, E{N), with 
three sets of Too bases constructed by stepwise inclusion of such elements in fully random 
sequences. We see that the convergence is not too fast; by a 50% truncation, we may lose 
200 keV of binding. 

Given the above results, it is interesting to see whether we can gain anything from a 
continuously random basis. While in the previous attempt the labels, which are natural 
numbers, were picked randomly, now the u values themselves, i.e. two real random num- 
bers were generated independently. In our procedure the values of both u were limited 
by imposing 

1 fm < z/"^/2 < 10 fm. (10) 

The bases generated in this way have proved excessively redundant. However, since the 
pool of basis states is infinitely large, it is reasonable to introduce some preselection. 
With this idea in mind, we introduced an admission procedure for extending the basis. 
Since any extension of the basis lowers the energy, the performance of a basis element 
can be judged by the energy improvement it causes. This obviously depends on the 
basis set in the previous step. At every step we tried 10 fully random basis states by 
adding one at one time in the basis constructed in the previous step, and admitted only 
the one that produced the lowest energy. A ((z^i^J"^^^, {t^2k2)~^^'^) obtained by such 
sequential random augmentations of the Too basis is dispalyed in fig. |^c, and the results 
of three such sets of calculations are shown in fig. |^. All three sets converge rapidly. At 
dimensions 30 and 50 the energies are —0.367, —0.380, —0.376 MeV and —0.393, —0.395, 
—0.394 MeV, respectively, while the corresponding radii are 2.403, 2.442, 2.458 fm and 
2.478, 2.467, 2.474 fm, respectively. Accepting an accuracy of, say 20-30 keV, we can say 
that, with such a random basis, a dimension of ~30-40 for a single subspace is adequate. 
The reduction of the dimension with respect to the discrete case is obviously due to the 
selective sampling. 

This method is comparable with stochastic basis selections that avoid full re-diagonal- 



ization at every step. Following ref. we repeated the above procedure by replacing 
the full diagonalizations with those of 2 x 2 matrices formed by the model g.s. obtained 
in the previous step and the state to be tested. As is seen in fig. ^, this procedure 
converges to the wrong value, while the same bases with full diagonalizations throughout 
yield a convergence almost as good as that resulting from the selection based on repeated 
full diagonalizations. Thus, whereas the limiting energy obtained by sequential 2x2 
diagonalizations is useless, the resulting basis is quite good. For problems in which the 
the diagonalization takes a higher proportion of the computing time, this procedure is 
recommendable as a method for basis construction, but one N x N diagonalization in the 
end is unavoidable. 

When the three subspaces are combined, the random selection with full re-diagonal- 
izations results in fig. Three of the four curves were calculated by choosing 40, 20 and 
20 basis elements from the subspaces consecutively, with three different sequences of the 
subspaces. It is seen that, while in the first subspace the convergence needs almost 40 
elements, in the others it suffices to take substantially fewer. By subtracting the lengths, 
along which the curves are practically horizontal, we see that the fastest convergence 
could be attained by starting with subspace Yn, which yields the deepest binding in a 
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single-subspace model (cf. table |l|). In calculating the forth curve, the channel of each 
basis element was also chosen randomly. The initial section of this convergence curve is 
substantially steeper than the others, but, as we see in fig. |^, this is not a gain if we 
prescribe an accuracy of 20-30 keV. With the random sequence of channels a 20 keV 
accuracy is produced by dimension 45 {E = —0.997 MeV). 

We have seen that a randomization with preselection may reduce the required dimen- 
sion appreciably, but, in general, a random selection does not produce the optimum. It 
is possible that one can get a basis that is closer to the optimum by including all those 
and only those random elements that surpass a certain preset utility level. Utility is not 
an abstract quality; it should be understood with respect to the basis generated in the 
previous step. This level can thus be defined as a minimum energy improvement e re- 
quired of the inclusion of a would-be basis state. When the admittance is decided by a 
utility test, the convergence is automatically signalled by insistent failure to find further 
elements that pass the test. 

In a pure Too model the utility test with e = 0.005 MeV did not improve the conver- 
gence, but in fig. lid we still show such a basis. This apparently differs from the preselected 
random basis in two respects: it does not contain highly overlapping pairs of basis states 
and it contains fewer elements with both small. In the full model this method looks 
definitely superior. In three runs, all with the sequence {Too,Yoo,Yii}, the energies and 
radii have converged to -0.999, -0.997, -0.998 MeV, and 2.429, 2.422 and 2.427 fm, at 
dimensions 45, 59, 47, respectively. When no basis element was found in a subspace to 
pass the utility test out of 10 consecutive attempts, the procedure was continued with 
the next subspace. In the last subspace, however, the calculations were only terminated 
when no basis element was found to fulfil the quality criterium after 100 attempts. The 
resulting energy and radius as a function of the dimension are shown in fig. |^. With 
£=0.01 MeV the energy deteriorated by a mere ~20 keV (it converged to —0.986, —0.984, 
—0.969 MeV at N =48, 39, 46, respectively), with a slightly larger scatter in the radii. 

Although the energy deepens monotonously and the rate of deepening is mostly found 
to decrease with the basis increased, there is no theorem on the latter. For this reason, 
due caution should be excercised to avoid declaring convergence prematurely. In view 
of this, we modified the method as follows. After A/" = 10 attempts failing to reach the 
utility level, we skipped the test once, i.e. we admitted a basis element that does not pass 



it. Such a procedure is called "simulated annealing" ||2^. The random search for new 
elements was stopped after 10A/'=100 consecutive failures to satisfy the utility criterium. 
With simulated annealing and parameter e = 0.01 MeV the energy curves converged 
to —0.993, —0.996, —0.992 MeV at dimensions 59, 57 and 66, respectively. Thus the 
simulated annealing makes the method less economical, but it does make the procedure 
that involves utility tests safer. 

Since the energy convergence attainable is not precise in a mathematical sense, one 
should check the quality of convergence for the wave function as well. The overall conver- 
gence can be judged from the fluctuations of the radius, which are acceptable for bases 
about to produce energy convergence (see figs. ^ and|^). The critical domain of the nu- 
cleus is, however, the halo region, which may be extensive without contributing too much 
to the binding. For example, /3-decay observables seem to be affected by the precise shape 
of the fall-off of the halo beyond 13 fm [23]. We tested the wave function through exam- 
ining, in the Too model, the two-variable function obtained by projecting out, from the 
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wave function (oo)o]oo' ^ intrinsic motion and the (/1/2) = (00) angular component 
of the relative motions: 

g[R^,R^) = R1R2 JdK,Y*,{K,) JdK2Y*o{K2){A{<^>oS{Ki - Pi)5(R2 - P2)}l^fo,(oo)o]oo)- 

(11) 

This is in fact a spectroscopic amplitude. We examined the effect of including, in the 
basis, elements that emphasize the surface region. Starting with the 100-dimensional basis 
reported on in table |l|, we included 10 more elements, resulting in a change of —0.00024 
MeV and +0.0027 fm in the energy and radius, respectively. Meanwhile, the norm square 
of g{Ri, R2) changed from s =1.3570 to 1.3568, and the change of g{Ri, R2) is very small 
accordingly. At the limit Ri = R2 = 15 fm, where g{Ri, R2) is 100 times smaller than its 
peak value, the change is 20%, and rapidly decreases with any of Ri smaller. Therefore, an 
integral over any of the two variables, like the function that enters into the description of 
f3 decay ^3], carries a much smaller error. Thus the error found seems insignificant even 



for the halo properties. A random basis constructed with utility tests (^=0.005 MeV) 
and simulated annealing converges at = 36, to yield an 0.008 MeV loss of binding 
and a —0.001 fm decrease of the radius (with respect to the 100-dimensional case) and 
s = 1.3596. The departure in g{Ri,R2) for large R values is, however, smaller than in 
the previous comparison, obviously because the difference between these two bases is not 
enhanced artificially in the asymptotic region. The slightly larger difference in s is due 
to a difference in g{Ri,R2) spread over a larger domain of (i?i,i?2)- Thus we can claim 
that the basis reduction does not hit the asymptotic region more than any other region, 
and the method is as reliable in predicting the halo properties as any other properties. 

To sum up, we have demonstrated that, by constructing the basis with stochastic 
sampling, the dimension required for an accuracy of 20-30 keV can be safely reduced by 
a factor of 5-6. We note that all results with random bases presented in this section have 
been corroborated by at least three different samplings. When there are more clusters, the 
basis is expected to be even more redundant, and an even larger reduction can probably 
be achieved. At the same time, however, there is more chance to be bogged down in flat 
areas, thus utility testing should always be accompanied with simulated annealing. 



4 Description of He 

Since in the a+n+n+n+n model considered there are four Jacobi coordinates, it is beyond 
our means to construct a basis with a prescription like eq. (^) even in a single arrangement. 
(Then the dimension would be 10^.) The experience with ®He suggests, however, that, if 
numerical convergence can be reached with a random basis, then there is good hope that 
the limiting energy and wave function are approximately equal to the exact energy and 
wave function of the model problem. Moreover, it is easy to ascertain that this is so by 
examining whether different random bases converge to the same energy. 

Fig. ^ shows that calculations with (preselected) random bases limited by eq. (|TD|) 
converge rather slowly. The cause of this is to be found in that the domain of the four- 
dimensional 'V-space" included according to eq. (|T^) is rather large in comparison with the 
domain that gives rise to net attraction between the five clusters. By squeezing the limits 
in eq. (|lOl), the convergence can only be speeded up for small dimensions because the outer 
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regions of this volume also contribute to the binding appreciably through interference in 
a configuration mixing. Therefore, we kept eq. ( [T0| ) as it is. 

Fortunately, however, the convergence can be speeded up by utility-tested basis states 
more radically than for ^He. In the example of fig. in subspace {[Q;(nn)]n}n with all 
orbital angular momenta zero, a preselected random basis leads to an apparent energy 
convergence at E = —2.970 MeV (measured from the a+n+n+n+n five-body breakup 
threshold) at dimension N = 300. Almost the same energy, E = —2.961 MeV, can be 
reached by utility testing of quality parameter e = 0.005 MeV, with a simulated annealing 
of A/" = 10, at dimension = 151. Therefore, we adopt this prescription as the standard 
procedure. The energy is accepted as converged when failure to pass the test is detected for 
the 100th time consecutively. In multiconfiguration calculations we limited the dimension 
for each subspace to 120, which was either enough to attain convergence in the above 
sense or was very close to convergence. In the latter case it was left to the configurations 
to be included afterwards to make up for this slight omission as far as this is within their 
scope. We checked all convergence in all cases by at least one more independent set of 
calculations, and adopted the lower energy as the calculated model energy. The agreement 
between such independent calculations has always been better than 30 keV in energy and 
0.005 fm in radius. 

A review of the results is given in table In all configurations considered the neutron 
spins are coupled to zero pairwise; in setting up these pairs we followed the pattern of the 
Jacobi vectors. All orbital angular momenta are assumed to be zero except those belonging 
to {[(an)n]n}n', which are coupled as {[(^i, ^2)^12, ^3]^i2,3, ^4}-^^ = {[(1, 1)0, 0]0, 0}0. The 
configurations considered are chosen so that each qualitatively different arrangement be 
represented. 

If all angular momenta were included in each arrangement, each arrangement would 
span a complete basis by itself. Therefore, a comparison of the single-configuration cal- 
culations can be interpreted as testing the performance of the single-angular-momentum 
ansatz in each arrangement. Alternatively, since each configuration may as well be viewed 
as exhibiting certain types of intercluster correlation, these calculations may be interpreted 
as exploring these correlations. The difference between the configurations is surprisingly 
small. The good energy of the shell-model-like configuration with pure / = shows the 
enormous flexibility of this model; after all, in the pure shell model a purely I = con- 
figuration would carry an excitation of Ahul The fact that the a-|-(4n)-like configuration 
Q;{[(nn)n]n} performs so moderately shows that, in the "all / = limit", the neutrons 
are more closely co-ordinated with the a cluster than with each other. At first sight, it 
looks surprising that the two Q;-|-(2n)-|-(2n)-like configurations do not work equally well; 
Q;[(nn)(nn)] is less satisfactory obviously because it is closer to an a-|-(4n)-type configura- 
tion like a{[(nn)n]n}. The arrangement (an)[(nn)n] is a representative of a-|-(3n)-|-n-like 
formations, and, not surprisingly, is energetically disfavoured with / = for similar rea- 
sons. The a;-|-(2n)-|-n-|-n-like arrangement yields a deep binding, which suggests that a 
configuration with two closely co-ordinated and two loosely co-ordinated neutrons, which 
is consistent with a two-neutron halo, has a large weight. The best result is, however, 
furnished by the {[(Q;n)n]n}n' configuration, which is not entirely unexpected since it 
contains, as it were, the best single ^He configuration, Yn. 

Some of the multiconfiguration calculations show that the single-configuration models 
leave some room for improvements. The component a{[(nn)n]n} does not seem to add very 
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much to the best single configuration with I = 0, but the combination of the two favoured 
configurations, {[a(nn)]n}n and [a(nn)](nn), do improve on the energy substantially. The 
third favoured / = configuration, {[(an)n]n}n, in {[a;(nn)]n}n+[a(nn)](nn)+{[(an)n]n}n 
still has an appreciable effect. Pilot calculations show, however, that a forth configuration 
has no significant effect. In particular, although the configuration containing / = 1, 
{[(an)n]n}n', is the best by itself, it scarcely deepens the binding when added to the 
superposition {[a(nn)]n}n+[a(nn)](nn)+{[(an)n]n}n. This statement lacks preciseness 
because the / 7^ ingredients make full-fiedged computations of such types prohibitively 
slow. We managed to include 17 basis states of subspace {[(a;n)n]n}n' and that improved 
the energy by a mere 0.0015 MeV. 

It was shown in ref. that it is possible to reproduce the energies of ^Li and ^Be 
in the cluster model with the same central force rather well, but ^He is much more 
complicated than ^Be. Now we see that, just as ^He, the nucleus ^He is somewhat 
overbound by the force chosen. By modifying the mixing parameter from u = 1.14 
to 1.135, the binding of ^He reduces to -0.940 MeV (with r = 2.462 fm), while the 
{[«(nn)]n}n+[a(nn)](nn)+{[(a;n)n]n}n model gives —3.041 MeV (with r = 2.272 fm) for 
^He. A linear interpolation shows that the force that puts the energy of ^He right yields 
—3.204 MeV for ^He, which is in excellent agreement with experiment, indeed. It is also 
remarkable that, while the separation energy of '^He changes significantly, there is virtu- 
ally no change in the radius. This suggests that the energy change, and thus the energy 
itself, is smeared all over the single-particle degrees of freedom. 



5 Conclusion 

We have argued that a thorough understanding of neutron-halo nuclei calls for a micro- 
scopic multicluster approach. In this paper we have demonstrated that such an approach 
is feasible. For the time being, we have adopted a central nucleon-nucleon interaction and 
a pure h.o. configuration for the a particle, and described the g.s. of the a+n+n and the 
a+n+n+n+n systems. We formulated the problem in a model space consisting of a few 
arrangements of the a and single-neutron clusters moving with orbital angular momenta 
of or 1 around each other. This assumption takes shape in well-defined model state 
spaces, and we wished to get nearly exact g.s. solutions in these spaces. Our approach 
uses a linear variational method, and the main technical problem was the construction of 
a basis that contains this "exact" solution. 

These multicluster problems are few-body problems aggravated by the complications 
coming from the internal structure of a. The exact solution of a nuclear five-body problem 
is beyond the present-day technical limits, and a most successful model-free approximation 
to such problems applies stochastic techniques in the variational calculus combined with 
similar (Monte Carlo) integration techniques (see e.g. |^). The method we have now 
introduced for multicluster problems only applies random sampling techniques at the level 
of the variational approach, viz. in constructing the trial function. This approach is still 
feasible because all matrix elements involved can be calculated analytically, even though 
they need substantial amount of computation after all. 

Drawing on former results, for ^He we adopted a three-component state space, and 
constructed a basis in a systematic way. Invariance of the results against amendments to 
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the basis has proven that the basis contains the "exact" g.s. to a good approximation. 
Having such an "exact" solution, we tested ideas of constructing smaller bases that still 
encompass the "exact" solution. It was found that step-by-step enlargements of the basis 
by qualified random elements can be used to construct smaller but almost as good bases. 
What qualifies a basis function is its contribution to the energy when added in the basis 
constructed in the previous step. We tried to select candidates for basis states both by 
comparing them with each other (i.e. we singled out the best among a few in every step) 
and with an absolute scale (i.e. we singled out all that contribute more than a preset 
value). Since the latter type of selection ("utility testing") might terminate itself too 
soon, we allowed temporary lapses in the rigour of admittance ("simulated annealing"). 

For ^He the same stochastic methods of constructing compact bases converged, and 
different methods or different random sets led to the same limiting energies. We thus 
have good reason to believe that these limiting energies are in fact the "exact" energies. 
While for ^He any of the methods works equally well, for ^He the utility testing proved 
definitely superior. To describe ®He, it seems to suffice to include a three-component 
state space, which is a superposition of an a+(2n)+n+n-type, an Q;+(2n) + (2n)-type and 
a shell-model-like arrangement. 

In spite of the simphcity of the interaction used, we managed to reproduce the g.s. 
energies of ^He and ^He simultaneously with an accuracy of 100 keV, which is a remarkable 
success. Nevertheless, the calculated matter radii of both nuclei undershoot the most 



recent experimental data. This result is consistent with that of Csoto for He, who 
used non-central terms along with the same central force, and extended the state space 
accordingly. These radii may be accounted for by the a-particle ingredient being too small; 
for the a wave function involved in the multicluster model has the same size parameter 
as producing 1.363 fm for the radius of the free a particle, which is ~ 8% smaller than 
the empirical value. This and other physical implications of the model will be discussed 
in a forthcoming paper. 
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Table 1: Energies and rms radii of ^He in different model spaces 



Subspace 


Too 


Yoo 


Yii 


{Too, Yoo} 


{Too,Yii} 


{Yoo,Yii} 


{Too,Yoo,Yii} 


Experiment 


E (MeV) 
r (fm) 


-0.395 
2.455 


-0.387 
2.441 


-0.435 
2.451 


-0.944 
2.411 


-1.015 
2.443 


-0.605 
2.469 


-1.016 
2.442 


-0.975 
2.4 ^) 



) 2.48±0.03, ref. |l|; 2.33±0.04, ref. [0. 



Table 2: Energies and rms radii of ^He in different model spaces 



Subspace ^) N E (MeV) ^) r (fm) =) 



{[(Q;n)n]n}n 


153 


-2.509 


2.341 


{[(an)n]n}n'; {[(1,1)0,0]0,0}0 


130 


-3.033 


2.322 


{[a(nn)]n}n 


161 


-2.961 


2.305 


[a(nn)](nn) 


163 


-2.676 


2.278 


a[(nn)(nn)] 


176 


-1.167 


2.188 


(an)[(nn)n] 


220 


unbound 




a{[(nn)n]n} 


157 


-1.368 


2.196 


{[a;(nn)]n}n+[a(nn)](nn) 


211 


-3.292 


2.306 


{ [a(nn)]n}n+a;{ [(nn)n]n} 


223 


-3.027 


2.225 


{[a(nn)]n}n+[a(nn)](nn)+{[(an)n]n}n 


303 


-3.395 


2.271 


{[a(nn)]n}n+[a(nn)](nn)+a{[(nn)n]n} 


311 


-3.321 


2.317 



^) The neutron spins are coupled to pairwise, and all orbital angular momenta are zero 
except those belonging to {[(an)n]n}n', which are coupled as {[(^i, ^2)^12, ^3]^i2,3, ^4}-^ = 
{[(1,1)0,0]0,0}0. 
^) Experiment: -3.112 MeV 
*^) Experiment: 2.52±0.03, ref. 



18|; 2.49±0.04, ref. |19 
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Figure 1: Schematic diagrams depicting all possible sets of Jacobi coordinates for ^He (a) 
and some possible sets for ^He (b). The a cluster and the neutrons are represented by 
large and small dots, respectively. 

Figure 2: The basis elements of the subspace Tqo represented by full or open circles in 
the (61,62) = («)"^/^ (^^2 )"^/^) sheet. The basis of (z/f^J-^/^ = g.T x 1.4*^' fm {h = 
1,...,10; i = 1,2) with the open circles omitted in the chessboard-like thinning (a); a 
random subset (full circles) of the previous basis (b); a random basis with preselection 
(c); a utility-tested basis (d). Note the logarithmic scale! 

Figure 3: Random paths of energy convergence in the Too subspace on the pool of discrete 
basis elements of (i^ffe.)"^/^ = 0.7 x lA^' fm {ki = 1, 10; ^ = 1, 2). 

Figure 4: Convergence in the Too subspace (horizontal line: exact value) with (preselected) 
random bases in energy (a) and in the rms point matter radius (b). 

Figure 5: Convergence in the Tqo subspace (horizontal line: exact value) with consecutive 
2x2 diagonalizations (short dashes), with full diagonalizations on the same bases (long 
dashes) and with a random basis constructed with full diagonalizations (full line). 

Figure 6: Energy convergence with (preselected) random bases in the {Too,Yoo,Yii} space. 
The four curves differ in the sequence of subspaces switched on. Sequence Too,Yoo,Yii: 
short-dashed curve; sequence Yoo,Yii,Too: medium-dashed curve; sequence Yii,Too,Yoo: 
long-dashed curve; basis elements randomized between subspaces: full curve; horizontal 
line: exact value. 

Figure 7: Convergence of energy (a) and rms radius (b) produced by three random 
sequences of bases constructed with utility tests of parameter e — 0.005 MeV in the 
{Too,Yoo,Yii} space. Horizontal line: exact value. 

Figure 8: Convergence of the ^He energy with a (preselected) random basis (dashed curve) 
and with a basis constructed with utility tests of £ = 0.005 MeV and jV = 10 (full curve) 
in the subspace {[Q;(nn)]n}n restricted by all orbital angular momenta set to zero. 
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